Particle-Hole Asymmetry and Brightening of Solitons in a Strongly Repulsive BEC 
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We study solitary wave propagation in the condensate of a system of hard-core bosons with nearest-neighbor 
interactions. For this strongly repulsive system, the evolution equation for the condensate order parameter of the 
system, obtained using spin coherent state averages is different from the usual Gross-Pitaevskii equation (GPE). 
The system is found to support two kinds of solitons when there is a particle-hole imbalance: a dark soliton that 
dies out as the velocity approaches the sound velocity, and a new type of soliton which brightens and persists 
all the way up to the sound velocity, transforming into a periodic wave train at supersonic speed. Analogous to 
the GPE soliton, the energy-momentum dispersion for both solitons is characterized by Lieb II modes. 



Strongly correlated quantum systems pose some of the most 
difficult challenges at the forefront of fundamental physics. 
Recent theoretical and experimental work in the field of Bose- 
Einstein condensates (BEC) of atomic gases center around 
unveiling new phenomena that could lead to the understand- 
ing of various complexities of these systems. The possi- 
bility of tuning inter-atomic interactions via Feshbach res- 
onances allows one to manipulate nature, and study realis- 
tic and tractable quantum many-body models. Among the 
various models, a system of impenetrable bosons known as 
the hard-core boson (HCB) gas, is a paradigm. It was ana- 
lyzed exactly in one dimension by Girardeau yjl. It has also 
been used to explore quantum phase transitions ||2|] and trans- 
port characteristics of bosonic and spin-polarized fermionic 
atoms llsl. 

In this paper, we study solitary wave propagation in a HCB 
system with nearest-neighbor (nn) interactions |0| to obtain 
deeper insight into beyond-GPE dynamics in quantum many- 
body systems. A soliton or solitary wave is a localized non- 
linear excitation that travels with a constant speed, retaining 
its shape. Non-dispersive solitonic energy transport has been 
observed in a BEcH. These are the well known dark soli- 
tons that travel with speeds less than that of sound. Such 
particle-like transport is an active field with particular em- 
phasis on unveiling many-body characteristics that cannot be 
described with the approximate description provided by the 
GPE. The existence of dark solitons have also been shown 
in a one-dimensional HCB gas using Fermi-Bose mapping 
Jstl, and in a generaUzed mean-field theory where the cubic 
nonlinearity of the GPE was replaced by a quintic term|0]. 
Solitons in similar quintic models were further investigated 
with a periodic potential|7p, and also in the presence of a 
dipole-dipole interaction. 181] Additionally, various numerical 
investigations^ have analyzed the quantum dynamics of dark 
solitons to study the effects of quantum fluctuations and quan- 
tum depletion in a Bose-Hubbard model. 

Our formulation, based on the equation for the BEC order 
parameter obtained using spin coherent state averages il ill dif- 
fers considerably from earlier studies. In addition to not being 
restricted to one-dimension, the evolution equation for the or- 
der parameter contains all powers of the condensate density. 



In this non-GPE type equation, which incorporates quantum 
fluctuations and depletion, both particles and holes emerge as 
equal partners in the transport. In this regard, our method- 
ology bears some parallel with recent work ll2ll on weakly- 
interacting atoms where a non-GPE description emerges due 
to quantum fluctuations. 

In contrast to the weakly repulsive condensate described by 
the GPE which supports only a dark soliton ! 1 311 . we show that 
the HCB condensate supports two distinct types of soUtons: 
a dark soliton whose amplitude vanishes as its propagation 
speed approaches the Bogoliubov speed of sound Vg (like the 
GPE soliton), and a new species of soliton that exhibits bright- 
ening and persists all the way up to Vg. Hence we call this a 
persistent soliton. For propagation speeds above Vg, it evolves 
into a periodic wave train. The existence of this novel type 
of soliton that brightens the condensate profile is tied to the 
particle-hole population imbalance, a key parameter for the 
cross-over to a non-GPE behavior. 

It is instructive to start our analysis with the extended lattice 
Bose-Hubbard model in d dimensions, 

H = — '^^[th^^hjj^a+Vnjnjj^a\+^^ Unj{nj — l)—{^—2t)n. 

(1) 

Here, bj and bj are the creation and annihilation operators for 
a boson at the lattice site j, rij is the number operator, a labels 
nearest-neighbor (nn) sites , t is the nn hopping parameter, U 
is the on-site repulsion strength, and fi is the chemical poten- 
tial. To soften the effect of strong onsite repulsion, we add an 
attractive nn interaction {V > 0), (although our results will be 
valid for repulsive V). Such a term may mimic certain aspects 
of the long range dipole-dipole interaction, that are subject 
of various recent investigations. Isl The on-site term 2tnj is 
added, so that the terms with t reduce to the kinetic energy 
in the continuum version of the many-body bosonic Hamilto- 
nian. 

The HCB limit [U oo) of the strongly repulsive Bose- 
Hubbard model corresponds to the constraint that two bosons 
cannot occupy the same site. This can be incorporated in 
the formulation by using field operators that anticommute at 
same site but commute at different sites, thus satisfying the 
same algebra as that of a spin-i system. The system can be 
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mapped to the following quantum XXZ Hamiltonian in a mag- 
netic fjeldlfiotl by identifying bj with the spin flip operator S^, 
along with 7ij = 



]{Ue-fi)S^. (2) 



3,a 



Here Ue = {t ~ V)d, where d is the spatial dimensionality. 

The dynamics of the HCB system (|2]i is described by the 
Heisenberg equation of motion. 



This operator equation can be transformed into an equation 
of motion for the condensate order parameter rjj = (5*^) 
using spin-coherent state averages jl 111 , a natural choice for 
describing the inherent coherence in the condensed phase 
of HCB. Parameterizing the local order parameter as T]j = 



: smf 



we obtain the condensate number density p'^ 



|77jp = i sin^ 6*^ and the particle number density pj = 
(uj) = {SJ S^) ~ sin^(i6'_,). Hence, the condensate and 
particle number densities are related by 



(4) 



P^ 



Therefore, in contrast to the GPE, we now have pj ^ pj. Fur- 
ther, the resulting formulation encodes fluctuations and deple 
tion, as seen from the relations {S^'Sj') — (S'^)(S'^) 
and {n]) - (n^f = 

The Hamiltonian in Eq. ([T]i is invariant, up to a change in 
the chemical potential, under a particle-hole transformation, 
where the hole operators are the hermitian conjugates of the 
boson operators, and the hole density is pj^l — pj. Thus 
the condensate density pj is the product of the particle density 
and the hole density (from Eq. with particles and holes 
playing equal roles in determining the condensate fraction. As 
we now show, the particle-hole imbalance variable Sj = (1 — 
2pj) = pj — Pj plays a key role in the dynamical evolution of 
the system. 

The equations of motion for Sj' and Sj lead to the follow- 
ing coupled equations for r/j and 6j : 

ihr]j = {Ue - p)vj -tSjY^ rjj+a + ^Vrjj ^ Sj+a, 

a a 

ihSj = 2t ^ (?7j rj*_^_^ - 77* 7]j+a)- (5) 

a 

We now consider the continuum approximation of the discrete 
equations (|5]), useful in the limit when the order parameter is a 
smoothly varying function with a length scale greater than the 
lattice spacing a. In the limit where the number of particles N 
and the number of lattice sites L tend to infinity, with the fill- 
ing factor N/L fixed, the system is described by the conden- 
sate order parameter ?7(r) = ^ sin6'(r) e"^'''', which is cou- 
pled to the particle-hole imbalance variable S{r) — cos 9{r): 





FIG. 1: (color online) Comparison of numerical (solid line: black) 
and analytical (dashed line: red) solitary wave solutions, for po = 
0.5. The upper and lower panels correspond to i) = 0, 7 = 1 and 
V — 0.87, 7 — 0.5, respectively. The plots with crosses (blue) show 
the corresponding condensate density p'^. 
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FIG. 2: (color online) Solitary wave for po ~ 0.25 for total particle 
density (a,b) and condensate density (c,d), for = (a,c) and v = 1 
(b,d). Solid (dotted) plots correspond to if ) solitons. 
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V-[(1-<52)V0], 



(7) 



where ta^ = h'^/{2m) and Ve ~ Va?. 

In the small-77 limit where we neglect terms involving 
|77p"77 {n > 1), and retain only linear terms involving the 
derivatives of both the order parameter equation (|6]l and the 
continuity equation (|7]l reduce to the GPE form, with an effec- 
tive mean field interaction equal to Up- In the above equations, 
if we set (1 — (5) = 2p and use Eq. (|4|i, we can show that for 
p{r) < i, the GPE limit is obtained when w p, i.e., the 
order parameter corresponds to the condensate of particles, 
while for p{r) > i, the GPE is satisfied by the order parame- 
ter for the condensate of holes, namely 77*, with p'^ « p'\ 

The Bogoliubov spectrum associated with the small am- 
plitude modes of the HCB system is similar to the GPE fll ill 
case and determines the Bogoliubov sound velocity Vs = 
(t/gjo'^/m)^/^. However, as we discuss below, for the propa- 
gation of nonlinear localized modes, true GPE-type transport 
emerges only near half-filling and the crossover to non-GPE 
dynamics is controlled by the particle-hole imbalance. 



FIG. 3: (color online) Energy vs canonical momentum for po =0.25 
for nn interaction V = Q. Branches (a) and (b) correspond to the 
persistent soliton and the dark soliton, respectively. Curve (c) gives 
the corresponding dispersion for the GPE soliton with p = p%. The 
dashed tangents on the two branches demonstrate linear dispersion, 
with identical slopes. 

We now investigate unidirectional soliton propagation in 
the HCB system, described by Eqs. ^ and dTji. To this end, 
we look for traveling wave solutions for p with velocity v, 
of the form p(£) = po + f{z), where z = {x — vt)/a and 
Po denotes a uniform background density. The corresponding 
background particle-hole imbalance is Sq = (1 — 2pq). The 
function f{z) is then found to satisfy the nonlinear equation 



1-C 



: /(/-'^o) 
Po 



= f- 



./(/ - '^0) 
Pi 



where 7^ = 1 — with v = v/vs- Further, z = (z and ( 



(8) 



is the Bogoliubov speed of sound Vg measured in units of the 
zero-point velocity h/ma. 

Equation ([D is satisfied by an elliptic integral. In the spe- 
cial case where the bracketed term on the left-hand side is 
approximated by unity, there are two soliton solutions of 
the formliil 



±^S^TWpI cosh (z/r) - ,5o 



(9) 



Here, F = 2(j is the width of the soliton. These two spe- 
cial solutions are emblematic of the general solution to Eq. 
(H). We have numerically integrated Eq. dS) to find f{z), and 
compared it with the analytical solution (|9]l. We find close 
agreement between the numerical and analytical solutions for 
all values of v and po- This happens because the neglected 
terms become irrelevant for small / and also when ^ = 0, 
the condition which determines almost all the contributions to 
the localized modes. Fig. 1 with po = 1/2 shows the case 
when the agreement is at its worst. Away from half-filling, 
the numerical and analytic solutions are very close. 

The solutions (|9]l are valid provided C/g > (or < > V), be- 
cause the functions are unbounded for negative C/g. We 
point out that the interaction V simply scales the width of 
the soliton, without altering its profile. This leads to an in- 
teresting possibility that the existence of the soliton and its 
overall profile may remain unaffected by certain long range 
interactions. Isl] 

The distinction between the twin solitons, which encode 
particle-hole duality as f^{z,po) = /^(z,Pq), can be elu- 
cidated by computing the momentum ~ m J J^dx. 
Here J = {h/m)pcV<p (see Eq. (|2|) is the condensate 
current, which is obtained using the equation lfl4ll d4}/dz = 
vma{p^ — Po) / hpc- We get 



= ±{mv) 



2C 



■ tan 



{SI + 2p07^)V2 ^ So 



(10) 



For small v, P+ corresponds to the momentum of a particle 
of positive mass while P~ corresponds to the momentum of 
a hole, as the effective mass is negative. This justifies associ- 
ating [respectively, /^] as the solitary wave with particles 
[holes]. Further, (P+ — P^)/mv = t^C/^/po^ which is sug- 
gestive of a conservation principle. 

Particle-hole imbalance appears to be a key factor in deter- 
mining the characteristics of the solitary waves. As seen in 
Fig. 1, when the number of particles is equal to the number 
of holes, the two solutions for the particle density p{z) are 
dark and anti-dark lfisll solitons ( bright solitons on a pedestal) 
which are mirror images of each other In the corresponding 
condensate density p'^{z), the solutions are indistinguishable, 
resulting in a single GPE-like soliton, that flattens out at the 
sound velocity. In fact, within the approximations discussed 
above, Eq. ^ expressed in terms of the condensate density 
Pc{z) is identical to a GPE with healing length ^ renormalized 



A/Vl — A^, where the microscopic parameter A — Vg raa/h to i^Vl — A^. In other words, when the background consists 
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of equal numbers of particles and holes, the soliton solution 
for the condensate density of the HCB is like the usual soliton 
of the GPE, but with a renormalized healing length. 

When the number of particles differs from the number of 
holes, the condensate soliton corresponding to the solution 
is dark and flattens out at the sound velocity. However, the 
condensate soliton corresponding to the solution /+ brightens 
the condensate profile, as shown in Fig. 2. As v increases, 
the spatial extent of the disturbance above the background 
increases, and the solitary wave p'^{z) becomes completely 
bright at the speed of sound. More important, this soliton does 
not flatten out but persists even at the sound velocity. The sur- 
vival or persistence of the soliton of this sector is evident from 
the limiting functional form lflill of /+ as u — > 1. We find 
p+ ^ po + + Cz^), where C = 4(^^(5o/pg and hence 

the width of the soliton becomes independent of its speed of 
propagation. It should be noted that as po 0, the solutions 
approach a definite limit with C equals to 8?7e in units of zero- 
point energy (ft/ of' jm. In other words, the anti-dark solitons 
become bright solitons, reminiscent of the bright solitons dis- 
cussed in systems with attractive dipole-dipole interaction.!^ 

Another distinctive feature of the persistent soliton is that 
it is transformed into a periodic wave train when its veloc- 
ity exceeds the sound velocity. Indeed, inspection of Eq. (|9]l 
shows that at supersonic speeds, when 7 becomes imaginary, 
the soliton is transformed to a spatially periodic wave that ex- 
ists for velocities < w < ^^/(l — (5q). 

Finally, we compute the energy-momentum relationship of 
the solitons by integrating the equation dE/dPc = u to ob- 
tain the canonical momentum Pc, where E is the energy 113 1. 
The three plots in Fig. 3 show the dispersion relations for the 
/+ and /" solitons, and the corresponding GPE sohton lll7ll . 



Each plot shows linear dispersion at one end of the spectrum 
(at the momentum corresponding to w = Vs), and saturation at 
the other end (at the momentum corresponding to v = 0). The 
dark soliton has a linear dispersion near E ^ 0, Pc ^ 0), with 
the slope given by Vs, and it saturates to zero slope near the 
maximum value of the momentum. (The GPE soliton, with 
p = /9q, has a similar behavior). The persistent soliton ex- 
hibits linear dispersion with the same slope Vg near the maxi- 
mum value of the momentum. Comparing this with the exact 
bosonic low-lying excitation spectrum for the HCB gas dis- 
cussed by Lieb [llSll, we note that both soliton branches with 
bounded momentum intervals, mimic the main characteristics 
of his type II excitation spectrum, being linear at one end and 
saturating at the other. However, it should be noted that in the 
present context, it is the lower branch that can be designated 
as a type II 'hole' state. The higher energy upper branch asso- 
ciated with the persistent soliton is in fact a type II 'particle' 
state. This is different from Lieb's classification, where 'par- 
ticle' states were always type I with unbounded momentum. 
In summary, we have explored solitary wave propagation in 



a degenerate HCB gas, and have found novel features that are 
not present in the conventional Gross-Pitaevskii equation, as 
well as some others that are. In a lattice system, this cross- 
over occurs at half-filling, which corresponds to equal num- 
bers of particles and holes in the HCB system. At exact half- 
filling, we obtain GPE-like dark solitary waves. Away from 
half-filling, we find both dark and anti-dark solitary waves 
including a persistent variety which propagates with a non- 
vanishing amplitude right up to the speed of the sound. The 
forms of these solutions at t = represent two initial types 
of disturbance profiles that could evolve into these solitons in 
realistic quasi-one-dimensional systems. These aspects could 
be further explored in numerical simulation in first-principle 
quantum many-body calculations and may find experimental 
realization in a highly anisotropic cigar shaped trap. We hope 
that our results will also provide further stimulus to the study 
of solitons in quantum many-body systems. 
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